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ABSTRACT 

Long-term variations are well-known in Be X-ray binaries, and are attributed to non- 
orbital changes in the structure of the Be circumstellar (equatorial) disc. However, the 
timescales involved are so long (tens of days, to years) that systematic studies have 
been very restricted. The ^8 year MACHO monitoring of the Large Magellanic Cloud 
(LMC) therefore presents an ideal opportunity to undertake such studies of Be X-ray 
systems that lie within the monitored fields. Here we present MACHO observations of 
two LMC Be X-ray transients, EXO 0531-66 and H 0544-665, the light curves of which 
show substantial (~0.5 mag) long-term variations. However, our analysis shows little 
evidence for any periodic phenomena in the light curves of either source. We find an 
upper limit for detection of a short (1-100 d) periodicity in the V- and i?-band light 
curves of EXO 0531-66 of 0.041 mag and 0.047 mag semi-amplitude, respectively. The 
upper limits for the V- and i?-band data of H 0544-665 are 0.054 mag and 0.075 mag 
semi-amplitude, respectively. Both EXO 0531-66 and H 0544-665 become redder as 
they brighten, possibly due to variations in the structure of the equatorial disc around 
the Be star. Spectra of both sources show Ha emission; for EXO 0531-66 we find the 
emission varies over time, thereby confirming its optical identification. 
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1 INTRODUCTION 

High mass X-ray binaries (HMXBs) generally fall into 
two subgroups, supergiant HMXBs and Be X-ray binaries 
(Be/XRBs). The companion star in a Be/XRB is a Be (or 
Oe) star, with a typical mass of ~10-20 Mq. A Be star 
is denned as a non-supergiant early-type star which has at 
some time shown Ha emission. This Balmer emission, along 
with a significant infrared excess, is believed to originate 
from circumstellar material which forms an equatorial disc 
around the Be star. The Be star significantly underfills its 
Roche lobe, but is thought to be rapidly rotating to explain 
the equatorial disc (Slettebak 1987). The orbital periods are 
long, ranging from ~15 days to years. 

The light curves of Be/XRBs, and isolated Be stars, 
display non-orbital long-term modulations which vary on 
timescales of tens of days to years. The modulations are 
thought to be due to a phase of matter ejection from the star, 
the Be phenomenon (Slettebak 1987). Due to the timescales 
involved it is difficult to study these long-term modulations 
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using normal observing programs. Hence, the observations 
obtained by the MACHO project provide a database from 
which to investigate such systems. 

In a previous paper (Alcock et al. 2001) we studied ~5 
years of MACHO observations of the Be X-ray transient 
(Be/XRT) A0538-66. The analysis revealed a long-term 
modulation of 421 d, together with the previously known 
16.6 d orbital period. We attributed this long-term period 
to the formation and depletion of the equatorial disc sur- 
rounding the Be star. 

These results motivated us to search for long and short 
term periodicities in other Be/XRTs. In this paper we 
present analysis of the MACHO light curves of two other 
Be/XRTs, EXO 0531-66 and H 0544-665. 



2 MACHO OBSERVATIONS 

The MACHO project monitored the LMC from 1992 July 
to 1999 December, with the primary aim of detecting grav- 
itational microlensing of constituent LMC stars by inter- 
vening dark matter. The observations were made using the 
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Figure 1. The MACHO light curves in blue (top panel) and red (bottom panel) filters of the proposed optical counterpart of EXO 0531- 
66, which is the northern component of a close double. The V and R magnitudes have been calculated using the absolute calibrations 
of the MACHO fields. The three arrows near day 100 in the top and bottom panels indicate the observations taken during the X-ray 
outburst in 1993 (Haberl et al. 1995), the time of peak X-ray luminosity is indicated by the middle arrow in each case. The arrow near 
day 1500 in the top and bottom panels indicates the time of the X-ray outburst observed by Burderi et al. (1998). 



1.27 m telescope at Mount Stromlo Observatory, Australia. 
A dichroic beamsplitter and filters provide simultaneous 
CCD photometry in two passbands, a 'red' band (~6300- 
7600 A) and a 'blue' band (~4500-6300 A). The latter filter 
is a broader version of the Johnson V passband (see Al- 
cock et al. 1995a, 1999 for further details). The 'blue' and 
'red' magnitudes were transformed to Johnson V and Kron- 
Cousins R respectively, using the absolute calibrations of the 
MACHO fields. 

The images were reduced with the standard MACHO 
photometry code SODOPHOT, based on point-spread function 
fitting and differential photometry relative to bright neigh- 
bouring stars. Further details of the instrumental set-up and 
data processing may be found in Alcock et al. (1995b, 1999), 
Marshall et al. (1994) and Stubbs et al. (1993), and we note 
that the MACHO database is now in the public domain 
(Allsman & Axelrod 2001). 



3 EXO 0531-66 

EXO 0531-66 was discovered in 1983 during observations of 
the LMC X-4 region with EXOSAT (Pakull et al. 1985). It 
was only seen for a month at the end of 1983 (Pietsch, Rosso 
& Dennerl 1989), and had not been seen in the EINSTEIN 
survey of the LMC (Long, Helfand & Grabelsky 1981). The 
source was detected again in 1985 with the SL2-XRT ex- 
periment on the Spacelab2 mission (Hanson et al. 1989). 
The lack of detection of the source in EXOSAT observations 
made between 1983 and 1985 indicates that EXO 0531-66 is 
a recurrent transient. 



The proposed optical counterpart (Haberl, Dennerl & 
Pietsch 1995) is a Be star, which is the northern component 
of a close double. During a 4 year observing program of 
EXO 0531-66 with ROSAT, a strong outburst was detected 
in 1993 which lasted more than two months (Haberl et al. 
1995). The intensity of this outburst was comparable to that 
seen by EXOSAT in 1983, and a 13.7 s pulsation period 
was detected (Dennerl, Haberl & Pietsch 1996; Burderi et 
al. 1998), indicating that EXO 0531-66 is a Be-neutron star 
binary. 

Haberl et al. (1995) suggested that the observed out- 
bursts were consistent with being caused by the periastron 
passage of the neutron star in a Be/XRB with an orbital 
period of 600-700 d. However, they could not rule out that 
the outbursts occur irregularly during a phase of increased 
matter ejection from the Be star. Dennerl et al. (1996) con- 
strained the orbital period to lie in the range of 4-70 d, 
assuming that the rate of period change of the pulse period 
was predominantly caused by Doppler shifts. They found an 
orbital solution with P or b = 25.4 d, e~0.1 and i~50°, for an 
assumed mass for the compact object of 1.4 Mg, and 15 Mq 
for the Be star companion. 

The identification of the optical counterpart remains 
uncertain as both the northern and southern stars in the 
close (~4") double are early type stars. Based solely on the 
X-ray error circle (9"), either star could be the optical coun- 
terpart. Stevens, Coe & Buckley (1999) obtained a spectrum 
of the northern component which confirmed the presence of 
Ha emission. A low resolution spectrum taken two nights 
later shows H/3 also in emission. No spectrum has been taken 
of the southern component. 
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Figure 2. PDM (top) and LS periodograms (bottom) for the l^-band data of EXO 0531-66. The frequency range and resolution is 
6.67xl0~ 4 -1.0 cycle d _1 and 1.25xl0 -4 cycle d — 1 , respectively. The dotted line is the 99% confidence level. 



We searched the MACHO database for sources within a 
9" radius of the X-ray position. We find that the only source 
which displays long-term variations typical of a Be star is 
the previously proposed optical counterpart, the northern 
component of the double. 



high energy light curve (15-60 keV). The V- and i?-band 
light curves show a slight brightening of the source at the 
time of the outburst reported by Burderi et al. (1998), this 
could be related to the X-ray activity. 



3.1 Light curve 

We show in Fig. ji] the 'blue' and 'red' photometry of 
the northern star proposed as the optical counterpart of 
EXO 0531-66, from 1993 January 19 to 1999 December 30. 
During three consecutive observations of EXO 0531-66, be- 
tween March and May 1993, Haberl et al. (1995) found the 
source to be X-ray active. The maximum observed count 
rate occurred on April 27 1993 and was found to correspond 
to 2.0xl0 36 ergs -1 in the ROSATb&nA (0.1-2.4 keV) for d = 
50 kpc. The times of the three observations taken during the 
extended outburst are marked on Fig. |lj Unfortunately, the 
MACHO project did not obtain observations of the source 
on these exact dates. However, observations were made near 
to these times and the light curve shows a slight brightening 
of the source which could indicate an optical response to the 
X-ray outburst. 

Burderi et al. (1998) observed EXO 0531-66 on March 
13-15 1997 with BeppoSAX (also marked on Fig. [j]). 
The estimated X-ray luminosity at the time of observa- 
tion, converted into the ROSAT band (0.1-2.4 keV), was 
~2.3xl0 36 erg s" 1 . This value is close to the luminosity of the 
outburst measured by Haberl et al. (1995). In a low state 
the source is found to have a luminosity of < 10 35 ergs _1 
(Haberl et al. 1995). Burderi et al. (1998) concluded that 
the source was in outburst at the time of their observation. 
During the ~2 day observation the light curve of EXO 0531- 
66 in the medium energy band (1.0-10.5 keV) was relatively 
flat, but with random fluctuations up to a factor of ~2 on 
timescales of 30 s. The low energy light curve (0.1-1.8 keV) 
displayed flaring episodes which were also observed in the 



3.2 Period Analysis 

Haberl et al. (1995) suggested that the orbital period for 
EXO 0531-66 was in the range of 600-700 d. However, Den- 
nerl et al. (1996) proposed a much shorter period in the 
range of 4-70 d, with a preferred value of 25.4 d. To inves- 
tigate variability on these timescales we performed a period 
search of the V- and 7?-band datasets over frequency range 
6.67xl0" 4 -1.0 cycle d~\ with a resolution of 1.25xl0" 4 cy- 
cle d~ . Before the temporal analysis was performed we de- 
trended the data by subtracting a linear fit. We searched for 
periodicities in the light curves using a Lomb-Scargle (LS) 
periodogram (Lomb 1976; Scargle 1982) and a phase dis- 
persion minimisation (PDM) periodogram (see Stellingwerf 
1978). 

We show in Fig. || the period searches of the V-band 
light curve of EXO 0531-66, the results for the _R-band data 
were similar. The largest peaks in the LS and PDM peri- 
odograms occur at ~1200 d and ~600 d. We folded both 
datasets on these periods. The folded light curves indicated 
that the ~600 d peak is a harmonic of the ~1200 d peak. We 
note that the value of 1200 d is close to half the length of the 
data. The light curves of EXO 0531-66 are dominated by the 
maxima at MJD 2449500 and MJD 2450750. Therefore, the 
1200 d peak is most likely due to the non-periodic variations 
characteristic of Be stars, rather than a true periodicity. 

We find no evidence for any short-term (1-100 d) pe- 
riodicites in the light curves of EXO 0531-66. The peaks at 
1 d in the LS and PDM periodograms are probably due 
to the sampling of the datasets. We find an upper limit of 
0.041 mag semi-amplitude for the presence of any short- 
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Figure 3. The optical light curves in blue (top panel) and red (bottom panel) filters of H 0544— 665 from MACHO project observations. 
The V and R magnitudes have been calculated using the absolute calibrations of the MACHO fields. 
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term periodicity in the V-band data, and 0.047 mag for the 
i?-band data. 



4 H 0544-665 

H 0544-665 was discovered with the HE A 0-1 scanning mod- 
ulation collimator by Johnston, Bradt & Doxsey (1979). The 
brightest star in the error region was found to have a mag- 
nitude of B~16 (star 1 in Figure 6 of Johnston et al. 1979), 
and was later found to be variable (Thorstensen & Charles 
1980). van der Klis et al. (1983) found that the star was 
variable on a timescale of weeks and determined a spectral 
type of BO-1 V. By plotting Fas a function of B- V, van der 
Klis et al. (1983) found a colour-magnitude correlation in 
which the star becomes redder as it brightens. This is typi- 
cal of Be stars whose variability is due to variations in the 
circumstellar disc. Coe et al. (1997) suggested that, due to 
the lack of IR emission from the candidate proposed by van 
der Klis et al. (1983), another star was the optical counter- 
part (Star 22 in Coe et al. 1997). However, the Be nature of 
star 1 of Johnston et al. (1979) was subsequently confirmed 
by Stevens et al. (1999) from spectra of the source in which 
double peaked Ha and weak H/3 emission is evident. They 
therefore concluded that this star is the optical counterpart 
of H 0544-665. 



4.1 Light curve 

Fig. ^ shows the MACHO project observations of H0544- 
665 taken during the period 1993 January 19 to 1997 Jan- 
uary 30. The magnitudes were calculated using the absolute 
calibrations of the MACHO fields (see § [] for details regard- 
ing the acquisition and reduction). 



4.2 Period analysis 

We detrended the V- and 7J-band data of H 0544-665 by 
subtracting a linear fit. We then searched for periodicities 
in the detrended light curves. The temporal analysis was 
performed over a frequency range of 1.43xl0 _3 -1.0 cycle 
d _1 , with a resolution of 1.25xl0 -4 cycle d — l . 

The LS and PDM periodograms for the i?-band data are 
shown in Fig. ^. The results for the F-band data were less 
significant. There is one LS peak which lies above the 99% 
confidence level, however, it is not significant at the 99.9% 
level. We conclude that there are no significant periodicities 
present in the light curves of H 0544-665. We determine an 
upper limit for detection of any short-term (1-100 d) period- 
icity in the V-band light curve of 0.054 mag semi-amplitude, 
and 0.075 mag semi-amplitude for the i?-band light curve of 
H 0544-665. 



5 COLOUR VARIATIONS 

We constructed V/(V-R) diagrams for EXO 0531-66 and 
H 0544-665 to investigate the changes in colour as the 
brightness of the sources vary. From Fig. ^ we can see that 
as EXO 0531-66 and H 0544-665 brighten their colour be- 
comes redder. This confirms the colour-magnitude correla- 
tion found for H 0544-665 by van der Klis et al. (1983). This 
indicates that the variations in the light curves are likely to 
be due to variations in the structure of the equatorial disc 
around the Be star. This is due to the disc being redder in 
B-V (i.e. cooler) than the Be star (Janot-Pacheco, Motch & 
Mouchet 1987). Thus, the formation of the equatorial disc 
will increase the optical brightness of the system by the ad- 
dition of red light, or it will make the system appear fainter 
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Figure 4. PDM (top) and LS periodograms (bottom) for the ij-band data of H 0544— 665, frequency range and resolution is 1.43xl0 — 3 — 
1.0 cycle d — 1 and 1.25xl0 — 4 cycle d — 1 , respectively. The dotted line is the 99% confidence level, the dashed line is the 99.9% confidence 
level. 



by masking the Be star, behaviour which is dependent on 
the inclination of the system. 



6 SPECTROSCOPY 

Spectra of the northern component of the double proposed 
as the optical counterpart of EXO 0531-66 were obtained 
on the nights of 1998 December 9 and 1998 December 14, 
using the 1.9 m telescope at SAAO. The detector was a 
SITe CCD attached to the Grating Spectrograph, using a 
1200 line/mm grating, giving a resolution of 1 A. Two 900 s 
exposure spectra of the source were obtained on both nights, 
together with Cu-Ar arc spectra for wavelength calibration 
and dome flats. Spectra of H 0544-665 (Star 1 of Johnston 
et al. 1979) were obtained on the nights of 1998 December 
11 and 1998 December 12, exactly as for EXO 0531-66. 

The spectra were reduced using iraf in a completely 
standard way. The one-dimensional spectra were extracted 
from the two-dimensional image using a FWHM aperture to 
ensure optimum signal-to-noise ratio. As no flux standards 
were observed, the spectra are presented in raw counts only. 



6.1 EXO 0531-66 

The summed spectra of EXO 0531-66 for each night are 
shown in Fig. |^ and confirm the presence of Ha emission 
as seen by Stevens et al. (1999). It is clear that the Ha 
emission profile is changing with time, as shown by the fac- 
tor ~2.3 increase in flux from the first night to the second. 
Balmer emission in the spectra of Be/XRBs is believed to 
originate in the circumstellar material disc around the Be 
star (Stevens et al. 1999). The variations in Ha we observe 
are most likely due to variations in this disc. 



6.2 H 0544-665 

Due to the poorer quality of the spectra of H 0544-665 we 
summed all four spectra in order to increase the signal-to- 
noise (Fig. . Ha emission is evident in the combined spec- 
trum, again confirming the results of Stevens et al. (1999). 



7 DISCUSSION 

To investigate the origins of the brightness variations that we 
observe in EXO 0531-66 and H 0544-665 it is instructive to 
consider another Be/XRT, A0535+26. This source has been 
extensively studied since its discovery in 1975 (Coe et al. 
1975; Rosenberg et al. 1975). The system contains a neutron 
star with a 104 s spin-period in an elliptical orbit of period 
~110.3 d around its 09.7 Hie companion (Priedhorsky & 
Terrell 1983; Nagase 1989). 

Clark et al. (1999; hereafter C99) analysed 15 years of 
optical photometry of A0535+26 and compared the variabil- 
ity found with that of the X-ray data taken with BATSE 
over the same period. As the neutron star accretes material 
from the circumstellar disc around the Be star, a study of 
the long-term variability of the envelope can help in the un- 
derstanding of the source's X-ray behaviour. A0535+26 was 
found to be highly variable in the optical over the 15 year 
period. 

C99 performed period searches of the long-term optical 
light curve of A0535+26. The analysis failed to find any 
evidence for modulation of the circumstellar envelope at the 
orbital period identified from the X-ray light curve. However, 
periods of ~1400, ~470 and ~103 d were found in the optical 
light curve, but their cause and whether they are coherent 
over time is unclear. C99 suggest that the changes in the 
optical light curve are due to the variability in the emission 
measure of an optically thin circumstellar envelope emitting 
via free-free and bound-free emission mechanisms. As there 
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Figure 5. V vs V-R plot for EXO0531— 66 (top panel) and H 0544-665 (bottom panel). Note that both sources get redder as they get 
brighter. 



was no evidence for optical modulation at the orbital period, 
C99 concluded that it is unlikely that the periastron passage 
of the neutron star effects the mass-loss rate from the optical 
companion. We note however that Hutchings (1984) found 
that the H/3 emission from A0535+26 was modulated on the 
111 d period. 

We find little evidence for any periodicities in the light 
curves of EXO 0531-66 and H 0544-665. The orbital peri- 
ods of Be/XRTs are usually determined from the recurrence 
period of outbursts, particularly in X-rays. We found a sta- 
ble period of 421 d for A0538-66 (see Alcock et al. 2001) 
which we suggest is related to the formation and depletion 
of the equatorial disc around the Be star, and is certainly 
not orbital in origin. However, the long-term modulations 
we observe in the optical light curves of EXO 0531-66 and 
H 0544-665 are most likely due to non-periodic variations 
characteristic of Be stars. If stable long-term variations are 
present in the light curves of EXO 0531-66 and H 0544-665, 
they are occurring on much longer timescales than presented 
here. 

If it is assumed that the orbital plane of the neutron 
star and the equatorial plane of the Be star are coincident, 
then when the neutron star passes through the dense mate- 
rial close to the Be star at periastron, an X-ray flare should 
be observed. C99 tested this model using their simultaneous 
optical and X-ray data of A0535+26. If the neutron star is 
accreting directly from the outflowing circumstellar disc ma- 



terial, the optical and X-ray light curves should be strongly 
correlated. 

C99 did not find any positive correlation between the 
optical and X-ray data, and concluded that this indicates 
that the neutron star is not accreting directly from the stel- 
lar wind. C99 found that the X-ray outbursts seemed to oc- 
cur after a period in which the optical light curve was seen to 
fade. This anti-correlation could represent a discrete episode 
of disc-loss which occurs radially, the interaction of the neu- 
tron star with this material then triggers the X-ray emission. 
They also observed a case in which an optical outburst was 
observed with no accompaning X-ray outburst. We find a 
slight brightening in the optical light curve of EXO 0531-66 
at times of recorded X-ray outburst (see Fig [ij Haberl et al. 
1995, Burderi et al. 1998). The lack of a strong correlation 
between the optical and X-ray data could be related to the 
low optical to X-ray luminosities observed for the sources. 

Janot-Pacheco et al. (1987) showed that changes in the 
emission from the circumstellar disc lead to a correlation 
between V and (B-V) for A0535+26. C99 confirmed that 
a correlation exists where the system becomes redder as it 
brightens, owing to a greater contribution from the cooler 
disc. We find the same correlation for the MACHO data 
of EXO 0531-66 and H 0544-665. However, another LMC 
Be/XRT we have studied, A0538-66, becomes redder as it 
fades (see Alcock et al. 2001). As the V/ (B-V) correlation 
is dependent on the inclination of the system, this suggests 
that we are viewing A0535+26, EXO 0531-66 and H0544- 
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Figure 6. Combined spectra of the Northern component of the double proposed as the optical counterpart of EXO 0531-66 for the nights 
of 1998 December 9 (lower) and 1998 December 14 (upper). The spectra have been smoothed and offset for clarity. The Ha emission line 
is marked, note the increase in emission from the first night to the second. 
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Figure 7. Combined, smoothed spectrum of H 0544-665. The Ho emission line is marked. 



665 at lower inclinations than A0538-66 (see McGowan & 
Charles 2002). 

Clark et al. (1998) investigated the long-term variability 
of optical spectroscopy of A0535+26. The profile and equiv- 
alent width of the Ha emission line were found to vary con- 
siderably over the 7 year period of study. It varied between 
a single-peaked and double-peaked or asymmetric structure, 
and the relative intensities of the symmetric and asymmet- 
ric profiles varied over time. The asymmetries were charac- 
terised by an additional feature in the blue or red shoul- 
der, which was sometimes displaced far enough from the 
rest wavelength to appear as a separate peak. Clark et al. 
(1998) conclude that the variability in the Ha line reflects 
the changes occurring in the circumstellar envelope, and sug- 



gest that the asymmetries observed in the profile are due to 
the envelope being asymmetric in geometry or density. 



8 CONCLUSIONS 

We can confirm the optical identification of EXO 0531-66 
with the northern component of the close double from the 
variability in the light curve, and the Ha emission that is 
observed in the spectra. The variability in the emission may 
be due to changes occurring in the circumstellar envelope as 
for A0535+26 (Clark et al. 1998). Our H 0544-665 spectrum 
also shows Ha emission, but more spectra are required in 
order to investigate its variability. 

Haberl et al. (1995) proposed an orbital period for 
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EXO 0531-66 of 600-700 d. Dennerl et al. (1996) proposed a 
much shorter period for the source in the range of 4-70 d. We 
find no variability on these timescales in the MACHO data. 
We find that the upper limits for a short-term (1-100 d) 
periodic modulation to be detected in the V- and TJ-band 
light curves of EXO 0531-66 are 0.041 mag and 0.047 mag, 
respectively, and 0.054 mag and 0.075 mag, respectively, for 
H 0544-665. 

The colour changes for A0535+26, EXO 0531-66 and 
H 0544-665 are consistent with the systems becoming redder 
as they brighten owing to the greater contribution from the 
cooler equatorial disc ( Janot-Pacheco et al. 1987) . This pro- 
vides evidence that we are observing the systems at a lower 
inclination than A0538-66 (see McGowan & Charles 2002). 
We detected the orbital period in the optical light curve of 
A0538-66 (see Alcock et al. 2001) which is also seen in X- 
rays. C99 find no evidence for modulation of the optical light 
curve of A0535+26 on the X-ray-determined orbital period. 
Hence, the system inclination may determine whether the 
orbital period of a Be system is detectable optically, and 
therefore we may only be able to detect those of EXO 0531- 
66 and H 0544-665 in X-rays. However, the detection of the 
orbital period in both optical and X-rays for A0538-66 may 
be connected to its short orbital period (16.6 d). If this is 
true, it suggests that the orbital periods for EXO 0531-66 
and H 0544-665 are both longer than ~20 d. 
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